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Most integral equations do not have a closed form solution. However, they can often be 
discretized and solved on a digital computer. 

Proof of the existence of the solution to an integral equation by discretization was first 
presented by Fredholm in 1903. 

In general, integral equations can be divided into two families: 

1. When the unknown is in the integral only, the integral equation is called of the first kind. 

2. When the unknown is both inside and outside the integral, the integral equation is called 
of the second kind. 

For electromagnetic applications, we can have both scalar and vector integral equations. 

1 Scalar integral equations 

Let us consider the situation of Fig. 1: two regions are defined in space, region 2 bounded by 
the closed surface S and region 1 being all the remaining space, bounded by S and Soo, in which 
sources are located. 
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Figure 1: Definition of the geometry of the problem. 



It is already known that we can write: 



(V 2 + ^)$i(f) = J(f), 
(V 2 + fc|) $ 2 (f ) = 0, 



(la) 
(lb) 
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Section 1. Scalar integral equations 



and similarly for the Green's functions 

(V 2 + kf) g\(f, r) = —5(f — r) in region 1, (2a) 
(V 2 + k\) g 2 (f, r) = -S(r - f) in region 2. (2b) 

Upon performing (Eq. (la)^i(f, r')— Eq. (2b)3>i(f, f')), we get: 

[c/i(f, r ) V 2 <I>i(f) - $i(f) V 2 5i(f, f')] = J(r)gi(f, r) + S(f - f')$i(f) . (3) 

Upon integrating Eq. (3) over the entire volume and using the identity V • (gV<£ — ^Vg) = 
gV 2 & - ®V 2 g, we get: 

/ dvV ■[g 1 (f,f')<f> 1 (f)-<f> 1 (f)Vg 1 (f,f')} = / dvJ(f) g(f, f') + $i(f') . (4) 

By Gauss theorem, we reduce the left-hand side integral to a surface integral. Also 

dvJ(f) gi (f,f') = -<S> inc (f'). (5) 
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We therefore obtain 

-/ dsn-[ S i(f,f / )V$i(f)-$i(f)V 5 i(f,f / )] = -$inc(r) + $i(r'), f ' € V . (6) 

By invoking the radiation condition, the integral over Soo vanishes, leaving (and exchanging 
f and f ' so that primed coordinates correspond to sources and unprimed ones to observation) : 



$i(f) = ^inc(r) - / ds'n- bi(f,f')V'*i(f / ) - $i(f')V^i(f, f')] f G Vi 



(7) 



For f G V2, the wave equation has no source and therefore the integration of the delta 
function yields a zero value. Performing the same steps for this second case, we get the generic 
relation: 



^inc(r) - / ds'n ■ [gi{f,f')V^{f') - $i(f')V 9l (f,f')} = I $lW r G Vl (8) 

JS f G Vo 



$i(r) r G Vi 
This is directly evocative of Huygens' principle: 



• For r € Vi: the total field 3>i(r) is the sum of the incident field plus the field due to the 
surface currents on the surface S. 

• For f G V2: the surface source on S produces a field that exactly opposes & mc , yielding 
the extinction theorem. 

Applying the same reasoning to region 2, we write (where there is no incident field): 

f G Vi 

(9) 

<J> 2 (r-) r G V 2 



! ds'n- [g 2 (f, r) V'$ 2 (f / ) - $2(r) Vg 2 (f, f')] 
Js 
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Note that the sign reversal is due to the definition of the normal vector h which has to point 
outward from the surface. Here, since we use the same n as before, we have to take it as being 
negative. 

Eqs. (8) and (9) have four independent unknowns: 

$i, <J> 2 ; ft ■ V3>i, n • V$2 , (10) 

which can be related by the boundary conditions. Here, g(f, r') and h ■ V^(f, f') are the kernel 
of the integral equation. 

2 Vector integral equation 

For the sake of completeness, we shall write the vector wave equation as well, although we will 
not use it here directly. 

Considering the same situation as above, we know that the fields have to satisfy: 



V x V x Ei(r) - uj 2 ei^iEi{r) =iu^iJ{r) , (11a) 

V x V x E 2 (f) - uj 2 e 2 ^2E 2 {f) =0 , (lib) 

and the Green's functions: 

VxVx5i(f, r ) - w 2 ei/xi^i(f, r) = 7S(f - r) , (12a) 
V x V x <5 2 (f,f') - uj 2 e2fi2^2(f,f') = 76(f - f 1 ) . (12b) 

By the same technique as before (Jydv [Eq. (11a) • G\{f,f') — Eq. (12a) • E\{f)]dv) 

/ dv Vx V xE 1 (r)-Gi(f,f')-E 1 (f)-V x VxGi(f,f') =iu^i dv J(f) -Gi(f, f') -E^f) . 
Jv L J Jv 

(13) 

The left-hand side can be transformed into a surface integral (left as exercise) and the right-hand 
side written in terms of incident field, yielding: 

E^f') = E mc {f')+ I dan-([Vx^i(f)]xSi(f,f / )+^i(f)xVxSi(f,f')j, f ' G Vy . (14) 

Js+Svc I J 

By reciprocity of the Green's functions, we can transform 

h- [V x ^i(f)] x Gi(f,f') = h x [V x ^i(f)] -Gi(f,f') 

= iujiiG\{f,f') ■ h x Hi(f) . (15) 

In addition, for an unbounded homogeneous medium (V x G(f,f') is reciprocal): 

h ■ Ei(f) x V x G\(f, f')=hx Ei(f) ■ V x Gi(f, r) 

= -[V x Gi(f,f')] -n x -Bi(f). (16) 
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Section 3. Problem with the internal resonance 



For Green's functions that satisfy the radiation condition, Eq. (14) becomes: 

Ei(f') = E inc (f') + J dsn- |i^iGi(f,r ) • h x #i(f) - [V x Gi(f,f')] • n x £i(f) j . (17) 

We perform the same steps for the other region to eventually obtain (and again interchanging r 
and f'): 



^inc(r) + J s ds ' |^iGi(f,f')-nx H^f') - [V x Gi(f,f')] - nx £i(f')j = < 
~y ds' ^^^(r, f') • n x i^r') - [V x W 2 {f,f')} ■ h x E 2 {f')^ = < 



Ei{f) f £ Vi 

o fev 2 

(18a) 

f G Vi 

E 2 (f) f £ V 2 
(18b) 

Together with the boundary conditions 

h x Hi(f) = h x ^(r) i (19a) 
n x ^i(f) = n x E 2 (r) , (19b) 

this system can be solved for h x E\{f) and n x H\{f). 

Note also that Eqs. (18) can be written in terms of electric and magnetic currents, and 
magnetic Green's functions: 

J eq (r') = hx H^f') ; -M eq {f') = hx E^f') , (20a) 

so that for example: 

E mc (f) + [ ds' \iijfiid ei (f,f') ■ J(f')+^ mi (r,f') ■ M(f')\ = \ El{r) r ^ Vl (21) 
Js I J r £ V 2 

3 Problem with the internal resonance 

A question arises: is the integral equation equivalent to Maxwell's equations? Or asked differ- 
ently, if we solve the integral equation and Maxwell's equations, do we get the same solution? 

The answer is actually "no", that is they are not always equivalent to each-other. The 
problem comes from spurious solutions at frequencies corresponding to the eigenfrequencies of 
the cavity enclosed by the surface S. This problem is generally referred to as the "internal 
resonance of the integral equation" . 

However, this lack of complete equivalence between the physical problem and its defining 
integral equation is rather minor and infrequent phenomenon, and is therefore often tolerated 
in practice. 
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4 Scattering by a rough surface 

Let us consider the 2D problem (for which we shall use a scalar integral equation) depicted in 
Fig. 2. 
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Figure 2: Rough surface S separating two media. 
The integral equation in scalar form is given by: 



$inc(r)+ / dsn- [<f>(f')Vg(f,f') - g(r,r )V$(f')] 
JS 



3>(r) r eVo 
f € Vi 



(22) 



4.1 Dirichlet boundary conditions: EFIE 

For a TE wave (E = Ey) and PEC surface, the boundary condition is 

<£>(f) = 0, for f e S . 
The integral equation becomes, for f G 5 , f' E S: 



$inc(r) - / ds' #(f,f') n • W(f') = 
Js 



3>(r) 




(23) 



(24) 



This equation, in which $ represent the electric field, is referred to as the electric field integral 
equation (EFIE). 

Note that as f gets closer to the surface, 5>(f) — > (from the boundary condition) so that 
we do not need to distinguish between approaching the surface from one side or the other. In 
fact, we can unify the equations and write: 



*inc(^) - / ds' g(f, r) n ■ V$(f') = 0, f £ S ,r £ S . 
Js 



(25) 



In addition, g(r,r') has an integrable singularity as r — > r' . Let us consider the surface 
depicted in Fig. 2, with z = fix): 



ds = Vdx 2 + dz 2 = dx^l + , 



(26) 



such that the integral equation becomes: 



*inc(f ) = J a dx'J 1 + (^j g(x, f(x),x', f(x')) (n • V$(f')) , at z> = f{x') , (27) 



G 



4.2 Neumann boundary conditions: MFIE 



where we can limit Ax to [-L/2, L/2\. 
By letting 



V 1+ (crQ (n-^(r'))\ z/=f{xf) = U(x f ), (28a) 

$ inc (x, f(x)) = b(x) , (28b) 
K(x,x')=g(xJ(x),x',f(x')), (28c) 



we can rewrite the integral equation as 

rL/2 
-L/2 

which has to be solved numerically. 

Before doing this, we shall write the integral equation with Neumann boundary conditions. 



rL/2 

/ dx'K(x,x')U(x') = b(x), (29) 

J -L/2 



4.2 Neumann boundary conditions: MFIE 

This corresponds to a TM case with H = Hy with a PEC surface. In this case, the boundary 
condition is 

n-V$(f) = 0, (30) 

and the integral equation becomes 

$inc(r)+ / da' $(f)h-Vg(f,f') = { (31) 
Js I f eV\ 

In this case, it makes a difference if we approach the surface from the top or from the bottom. 
In fact 

$inc(r+) + / ds'$(f') h ■ Vg{f+, f') = $(f+) , (32a) 
Js 

3>mc(r-)+ / ds'$(r) h-Vg(f^,r) =0. (32b) 
Js 

These two equations seems inconsistent with one another. We shall show that in fact, they 
are actually consistent with each-other, due to the singularity of the Green's functions. 

We shall examine what happens when we let f approach the surface. Fig. 3 is an illustration 
of the situation at the immediate vicinity of the surface. The integral part of the equation can 
be written as: 

f ds'$(f') h ■ Vg(f, f') = PV [ ds'®(f')h-Vg(f,f')+ [ ds'®(f')h-Vg(f,f'), (33) 

Js Js Jpiece 

where PV denotes the principal value and 'piece' refers to the integration over the local domain 
shown in Fig. 3. For this integral, we use the local coordinates: 

ds' = dX', h = Z', etc (34) 
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a x , 

Figure 3: Zoom on the rough surface. 



The integral becomes: 

f ds'$(f')h-Vg(f,f') = lim lim 3>(f) f dX'J^—g(f,r) 

./piece J-a 



i piece 

Over the small piece, we have: 



\f-f'\ = y/X' 2 + {Z- Z'f , 



g =-H { i\\f-f'\) - -£ln( 



r — r 



2vr 



<9# 



8Z> J z , =0 2tt(X' 2 + Z 2 )' 
Thus, the integral becomes 

f =lim lim $(f) / dX'- 

Jpicce 1^1-0 7_ a 27 

$(f) 



lim lim 

a^O |Z|->0 27T 



tan" 



2vr X' 2 + Z 2 



< 



I $(f) for Z > 
-2*(f) forZ<0 

The two parts of the integral then become: 



^inc(^) + PV [ ds'®(r ) h ■ Vg(f, r) + -$(f) = $(f) , 
Js 2 

^ine(r) + / ds'$(f') n • V 5 (f, f') - U(f) = , 

which is cast into one equation: 

^inc(r) +PV I ds'$(r) h ■ X7g{f, r) = ^$(f) . 
Js 2 

Eq. (39) is called the magnetic field integral equation (since $ represents the maj 
and is an integral equation of the second kind. 



(35) 

(36a) 
(36b) 

(36c) 



(37) 



(38a) 

(38b) 
(38c) 

(39) 

metic field) 
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Section 5. Solving the EFIE 



In the same way as before, we can define the kernel as 



K(x,x') = Jl+(^j h-Vg(r,r')\ z=f{xlz , =nx , h (40) 



and write the MFIE as 



roc 1 

$inc(f) + PV dx'^(x')K(x,x') = (41) 

J-oo 2 

5 Solving the EFIE 

Upon using the notation introduced in Section 4.1, the problem is reduced to solving the following 
integral equation: 

,L/2 

/ dx'K(x,x')U(x') = b(x). (42) 

J-L/2 

Let us subdivide the integration domain into N small elements, each of length A = L/N, 
and centered at x rn (me [1,N]). Thus, constraining the observation at these discrete locations, 
the integral equation becomes 



f L/2 

/ dx' K(x m , x') U(x') = b(x m ) . (43) 

J-L/2 



Next, if we suppose that U(x) is constant in each interval, we replace the integral by a sum over 
all segments, excluding the singular term: 

Ax y K(x m , x n )U(x n ) + U(x m ) dx'K(x m ,x') = b{x m ). (44) 

n=l Jm 

Note that we have to single out the singularity, i.e. the m th interval because K(x m ,x') is 
singular at x' = x m . This part is known as the self-patch contribution. 
For a 2D problem, we have: 

K(x m , x') = - A H^(ky/(x m - x'Y + (f(x m ) - /(x')) 2 ) . (45) 
Upon approximating f(x') — f(x m ) ~ f'(x m )(x' — x m ), we write: 



K(x m ,x') = % -H<i\k\x' -x m yi + f'{x m y). (46) 



For small argument: Hq (a) ~ 1 + z-ln(a7/2) where 7 is the Euler constant (7 ~ 1.78). 
Therefore: 

i 



1 + i X - lnQfc(x' - x m ) v 7 ! + / ,2 (^ m )) 



(47) 
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Therefore, the integral becomes: 



p rx m +Ax/2 rAx/2 

/ dx'K(x m ,x') = 2 / dx'K(x m ,x') = 2 / dx'K(x m , x + x m ) 

J m J x m J 



~ - / arc 
2 7o 



1 + i-ln( JfexVl + / /2 (^r 

7T V 2 



4 [ 7T \4e 
where ln(e) = 1. 

We can therefore cast the integral equation into a matrix equation of the form: 

N 



where 



^ ^ A m n U n — b m j 



n=l 



(48) 



(49) 



U n = U (x n ) is the unknown, 



Ax if(x m ,x r 
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1 + zfln [t Axy/1 + 



for m ^ n 
for m = n 



(50a) 
(50b) 

(50c) 



The system can now easily be solved numerically 



